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We investigate how field-effect doping affects the structural properties, the electronic structure 
and the Hall coefficient of few-layers transition metal dichalcogenides by using density-functional 
theory. We consider mono-, bi-, and trilayers of the H polytype of M 0 S 2 , MoSe 2 , MoTe 2 , WS 2 , 
and WSe 2 and provide a full database of electronic structures and Hall coefficients for hole and 
electron doping. We find that, for both electron and hole doping, the electronic structure depends 
on the number of layers and cannot be described by a rigid band shift. Furthermore, it is important 
to relax the structure under the asymmetric electric field. Interestingly, while the width of the 
conducting channel depends on the doping, the number of occupied bands at each given k point 
is almost uncorrelated with the thickness of the doping-charge distribution. Finally, we calculate 
within the constant-relaxation-time approximation the electrical conductivity and the inverse Hall 
coefficient. We demonstrate that in some cases the charge determined by Hall-effect measurements 
can deviate from the real charge by up to 50%. For hole-doped MoTe 2 the Hall charge has even the 
wrong polarity at low temperature. We provide the mapping between the doping charge and the 
Hall coefficient. In the appendix we present more than 250 band structures for all doping levels of 
the transition-metal dichalcogenides considered within this work. 


I. INTRODUCTION 

Since the rise of graphene^ and the discovery of 
topological insulators^ a lot of interesting physics has 
been found in systems with reduced dimensions. Other 
two-dimensional (2D) material, such as monolayers 
or few-layer systems (nanolayers) of transition-metal 
dichalcogenides^— (TMDs) are gaining importance be¬ 
cause of their intrinsic band gap. TMDs are MX 2 -type 
compounds where M is a transition metal ( e.g ., M = 
Mo, W) and X represents a chalcogen (S, Se, Te). These 
materials form layered structures in which the different 
X-M-X layers are held together by weak van der Waals 
forces. Thus, similar to graphene, one can easily extract 
single or few layers from the bulk compound using the 
mechanical-exfoliation or other experimental techniques. 

Doping these nanolayers with field-effect transistors 
(FETs) is particularly appealing^— as it allows for the 
exploration the semiconducting, metallic, superconduct¬ 
ing, and charge-density-wave regimes in reduced dimen¬ 
sionality. Furthermore, the TMDs are promising mate¬ 
rials to realize valleytronics, i.e ., the usage of the valley 
index of carriers to process information^ ' 15 ' 22 . 

Despite these challenging experimental perspectives, 
the understanding of structural, electronic, and transport 
properties at high electric field in FET configuration is 
still limited, particularly in the physically relevant case of 
multilayer samples. Previous theoretical work s 20 ' 23 ana¬ 
lyzed the high doping limit of 20 nm thick M 0 S 2 flakes^, 
relevant for ionic-liquid based FETs, by assuming that 
only the topmost layer is doped uniformly. However, it 
is unclear to what extent the doping of thick flakes can 
be modeled in this approximation as the thickness of the 
conductive channel is not experimentally accessible. 


A second crucial issue is the determination of the 
doping charge. Usually the charge is determined via 
Hall-effect measurements. However, the interpretation of 
Hall experiments assumes a 2D clectron-gas model, most 
likely valid only in the low doping regime. In TMDs, due 
to the multivalley electronic structure, this assumption 
is highly questionable. 

In this paper we solve these issues and provide a 
thorough study of structural, electronic and transport 
properties and of their changes under field-effect dop¬ 
ing for TMDs. We use our recently developed first- 
principles theoretical approach to model doping in field- 
effect devices^. The method allows for calculation of the 
electronic structure as well as complete structural relax¬ 
ation in field-effect configuration using density-functional 
theory (DFT). We apply our approach to the H polytype 
of M 0 S 2 , MoSe 2 , MoTe 2 , WS 2 , and WSe 2 - 

The paper is organized as follows: In section[TT]we sum¬ 
marize the parameters and methods used within the pa¬ 
per and the relaxed geometries of the bulk TMDs. Then 
we will first show the results for the undoped mono-, bi-, 
and trilayer systems (section IIII Al) . After a brief discus¬ 
sion on the quantum capacitance (IIII BI) . we will investi¬ 
gate the changes of the geometry and the electronic struc¬ 
ture of the different TMDs in sections IIII Cl and IIII Dl 
respectively. Finally, we will focus on how the doping- 
charge concentration can be determined experimentally 
by Hall-effect measurements. We will show in section 
IIII El that the results of such a measurement cannot be 
interpreted within a 2D electron-gas model but that the 
specific band structure of the TMDs and its changes in 
a field-effect setup need to be taken into account. In 
section IIII FI we will furthermore investigate the density 
of states (DOS) at the Fermi energy and the electrical 
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conductivity as function of doping. In the end, we will 
summarize our results and draw some final conclusion in 
section ED In the appendix we provide a full database 
of electronic structures for all doping levels considered in 
this work (in total more than 250 calculations, Figs. 28- 
63) and Hall coefficients (Figs. 64-66) of mono- , bi-, and 
trilayer dichalcogenides as a function of doping. 


II. COMPUTATIONAL DETAILS 

All calculations were performed within the framework 
of DFT using the Quantum ESPRESSO package^ 
which uses a plane-wave basis set to describe the valence- 
electron wave function and charge density. We employed 
full-relativistic, projector augmented wave potentials— 
While the local-density approximation is known to un¬ 
derestimate the lattice parameters, generalized-gradient 
approximations for the exchange-correlation energy over¬ 
estimate the out-of-plane lattice constant (see, e.g., 
Ref. H3 and references therein). In our FET setup a 
correct description of the interlayer distance is however 
very important. Accordingly, we choose the Perdew- 
Burke-Ernzerhof functional^ (PBE) for the exchange- 
correlation energy and furthermore included dispersion 
corrections^ (D2). This also leads to the best agree¬ 
ment with both the experimental in-plane and out-of- 
plane lattice parameters (c/., Tab. Q]). A comparison be¬ 
tween PBE and LDA for M 0 S 2 can be found in the ap¬ 
pendix (Fig. 27). 

Using the experimental lattice parameters of the H 
polytype of the bulk structures^ as starting geometry, 
we minimized the total energy as a function of the lat¬ 
tice parameters until it changed by less then 2 meV. For 
the molybdenum-containing dichalcogenides we used a 
cutoff of 50 Ry and 410 Ry (1 Ry ft 13.6 eV) for the wave 
functions and the charge density, respectively, while for 
the tungsten dichalcogenides we chose 55Ry/410Ry. The 
Brillouin zone (BZ) integration has been performed with 
a Monkhorst-Pack grid^i of 16 x 16 x 4 k points and using 
a Gaussian broadening of 0.002 Ry ft! 27 meV. The con¬ 
vergence with respect to the number of k points as well 
as the wave-function and charge-density cutoff has been 
checked. The self-consistent solution of the Kohn-Sham 
equations was obtained when the total energy changed by 
less than 10 -9 Ry and the maximum force on all atoms 
was less than 5 • 10~ 4 Ry ag 1 (do ft 0.529177 A is the 
Bohr radius). The lattice parameters thus determined 
are given in Tab. [I] and agree within 2% with the experi¬ 
mental values. 

The final geometry of the bulk system was used as the 
starting geometry for the relaxation of the layered-2D 
systems in an FET setup. To achieve this, the size of 
the unit cell was fixed in plane and the perpendicular 
size was increased such that the vacuum region between 
the repeated images was at least 23 A. The layers were 
stacked as in the H polytype of the bulk compound. The 
BZ integration has been performed with a Monkhorst- 


TABLE I. Comparison of the calculated and experimental 
lattice parameters. 



a calc. 

a exp. 

c calc. 

c exp. 

M 0 S 2 

3.197 A 

3.160 A 

12.38 A 

12.29 A 

MoSe2 

3.328 A 

3.289 A 

13.07 A 

12.93 A 

MoTe2 

3.536 A 

3.518 A 

14.00 A 

13.97 A 

WS 2 

3.190 A 

3.153 A 

12.15 A 

12.32 A 

WSe 2 

3.341 A 

3.282 A 

12.87 A 

12.96 A 


Pack grid of 64 x 64 x 1 k points for the charged systems 
and 16 x 16 x 1 k points for the neutral ones. In order 
to correctly determine the Fermi energy in the charged 
systems, we performed a non-self-consistent calculation 
on a denser k-point grid of at least 90 x 90 x 1 points 
starting from the converged charge density. All other 
parameters were the same as in the calculations for the 
bulk systems. For the total-energy calculations of the IT 
and IT’ polytype of M 0 S 2 and WSe 2 under FET doping, 
we doubled the unit-cell size along one in-plane direction 
and correspondingly halved the number of k points along 
this direction. 

Several methods have been developed to study electro¬ 
statics in periodic boundary condition™™ - — with differ¬ 
ent experimental geometries. We used our recently de¬ 
veloped method^! as it is tailored for the FET setup and 
allows for structural optimization. The dipole for the 
dipole correctio n 24 ! 35 was placed at Zdi P = ddi P /2 with 
ddip = 0.01 L and L being the unit-cell size in the direc¬ 
tion perpendicular to the 2D plane - L changed for the 
different calculations and was between 34A and 48A. The 
charged plane modeling the gate electrode^ was placed 
close to the dipole at z mono = 0.011L. A potential barrier 
with a height of Vq = 2.5 Ry and a width of <4 = 0.1 L 
was used in order to prevent the ions from moving too 
close to the gate electrode. The final results were found 
to be independent of the separation of the dipole planes, 
as well as the barrier height and width as long as it is 
high or thick enough to ensure that the electron den¬ 
sity at the position of the dipole and the gate electrode 
is zero, p e (^ mon o) = p e (zdi P ) = 0. As we will often 
give the doping-charge concentration per unit cell n (i.e., 
in charge per unit cell, e/unit cell, with the elementary 
charge e ft 1.602 x 10 -19 C), we summarized in Tab. El 
the conversion to charge-carrier concentration per area n 
(in cm -2 ) for the different dichalcogenides and two typ¬ 
ical doping-charge concentrations of n = 0.01 e/unit cell 
and n = 0.15 e/unit cell. Throughout the paper we will 
use n < 0 and n > 0 for electron and hole doping, re¬ 
spectively. Typical charge-carrier concentrations which 
can be achieved in experiments using either solid-state 
dielectrics such as SiC >2 or ionic-liquid based FETs can 
be found in Tab. m Note that in principle those max¬ 
imum concentrations could be possible for all TMDs^ 
even if we did not find references for, e.g., ionic-liquid- 
based field-effect doping of MoSe 2 - 
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TABLE II. Conversion for the doping-charge concentration 
n (in e/unit cell) to charge-carrier concentration n per area 
cm -2 for the different dichalcogenides and two typical doping- 
charge concentrations of n = 0.01 e/unit cell and n = 
0.15 e/unit cell. 



n 

= 0.01 e/unit cell 

n 

= 0.15 e/unit cell 

M 0 S 2 

n ft 

0.1127 • 

10 14 cm -2 

n ft 

1.6911 • 

10 14 cm -2 

MoSe 2 

n ft 

0.1042 ■ 

10 14 cm" 2 

n ft 

1.5636 • 

10 14 cm" 2 

MoTe 2 

n ft 

0.0924 • 

10 14 cm" 2 

n ft 

1.3856 • 

10 14 cm" 2 

WS 2 

n ft 

0.1135 • 

10 14 cm -2 

n ft 

1.7027 • 

10 14 cm" 2 

WSe 2 

n ft 

0.1035 ■ 

10 14 cm" 2 

n ft 

1.5521 • 

10 14 cm" 2 


TABLE III. Maximum experimental charge-carrier concentra¬ 
tion n per area cm -2 for the different dichalcogenides using 
either a solid-state or an ionic-liquid based FET. References 
are given after the charge-carrier concentration. 



Polarity 

Solid-state FET 

Ionic-liquid FET 

MoS 2 

n-type 

« -3.6- 10 13 , [14] 

ss -9.8 • 

10 14 , [37] 


p-type 

— 

— 


MoSe 2 

n-type 

« -3.9- 10 12 , [38] 

— 



p-type 

« +2.0- 10 12 , [38] 

— 


MoTe 2 

n-type 

« -1.3- 10 13 , [39] 

ss - 1.1 • 

10 13 , [40] 


p-type 

— 

«+1.8- 

10 13 , [40] 

WS 2 

n-type 

« -1.0 ■ 10 14 , [19] 

« -4.0 • 

10 14 , [41] 


p-type 

— 

« +3.5 ■ 

10 13 , [42] 

WSe 2 

n-type 

— 

« -1.4 • 

10 14 , yj 


p-type 

« +9.0 • 10 12 , [43] 

«+1.9- 

10 15 , [16] 


In order to calculate the Hall tensor Rijk(T; Ef), we 
used the BoltzTraP code!! to determine the conductivity 
tensors in Eqs. m and (fl2l) (see section UlI El for more 
details). We fitted the band structure for each doping of 
the different TMDs by using 55-times more plane waves 
than bands and used afterwards the in-plane components 
of the energy-projected tensors to calculate the Hall coef¬ 
ficient R X yz(T\ Ep) for temperature T and chemical po¬ 
tential Ep. R xyz {T\Ep) is the only relevant Hall coef¬ 
ficient for our 2D systems assuming that the magnetic 
field is applied perpendicular to the layers. We checked 
the convergence by calculating R xyz (T ; Ep) with increas¬ 
ing number of k points and found that the results for 
the 64 x 64 x 1 grid and the dense grid of the non-self- 
consistent calculation are the same. 


III. RESULTS 

A. Electronic structure of TMDs 

In the following we will first briefly summarize 
the results for the undoped TMDs before investigat¬ 
ing the changes under field-effect doping. We fo¬ 


cus on the changes in the valence-band maximum and 
the conduction-band minimum with changing transition 
metal or chalcogen and compare them with other results 
found in literature!!. 

Figure Q] shows the band structure and the projected 
density of states (pDOS) for monolayer MoSe 2 with and 
without including spin-orbit coupling (SOC). Monolayer 
molybdenum diselenide is (as most TMDs) a direct-band- 
gap semiconductor with a DFT gap of about 1.329 eV at 
the K point. Our calculated gap is smaller by 83 meV 
than the one in Ref. [2j] which can be attributed to our 
slightly larger in-plane lattice parameter as the size of the 
band gap decreases with increasing lattice constant!! - — . 
This is due to the fact that the valence-band maximum 
at K is formed by in-plane states of both the transi¬ 
tion metal and the chalcogen!!. On the other hand, the 
conduction-band minimum at K is mainly formed by out- 
of-plane Mo states (d z 2 without SOC and rrij = ± 1/2 for 
both j = 5/2 and j = ?>/2 including SOC) and in-plane 
states of the chalcogen. The valence-band maximum near 
T has basically only out-of-plane states of Mo and Se as 
can be seen in Fig. [T] This will become very important 
for hole doping of the nanolayers in an FET setup - de¬ 
pending on which valley is doped (K or T ) one can expect 
different doping-charge distributions. Energetically very 
close to the conduction-band minimum at K is a min¬ 
imum half-way between K and F. The corresponding 
point in k space is called Q in literature (sometimes A or 
A m in as it is a minimum along the A line, from T to K) 
even if it is not a high symmetry point of the BZ. This is 
also why this minimum does not lie exactly at the same 
point for the different TMDs and its position can even 
change if the number of layers is increased. The states 
close to the Q point have a stronger in-plane character 
and can thus also lead to a different doping-charge dis¬ 
tribution if the doping occurs mainly at this point in the 
BZ. The same results for the character of the different 
valleys were also obtained in Refs. l47l - |49i 

The different character of the states in the different val¬ 
leys is even more important for the other TMDs. From 
sulfur to tellurium the difference between the minimum 
at K and Q decreases: for M 0 S 2 the minimum at K is 
lower by 279 meV while it is only 154 meV and 72 meV 
lower for MoSe 2 and MoTe 2 , respectively. The change in 
the case of the tungsten dichalcogenides is much lower 
which however might be due to the stronger spin-orbit 
splitting of the bands near Q compared the splitting at 
K. The band structures of all undoped TMDs are sum¬ 
marized in the appendix, Figs. 19-23. 

Increasing the number of layers in TMDs leads to a 
well-known change from a direct-band-gap semiconduc¬ 
tor to an indirect on o 27 ' 49 i 51 — as shown in Fig. [2j The 
change of the direct band gap at K with increasing the 
number of layers is much smaller than the changes at 
r or Q. This is due to the small hybridization between 
different layers at K as those states have only in-plane 
chalcogen character. On the other hand, both valleys 
at r and Q have contributions from Se p z states. Ac- 


















4 


MoSe 2 


> 

<D 


k( 


> 

<u 

kf 

kJ 



_i_ !L __ W / ___I_,_i_ 

r m k r i 2 

DOS [states / eV] 


FIG. 1. (color online) Band structure and density of states 
projected onto atomic orbitals for monolayer MoSe 2 with¬ 
out (upper panel) and with (lower panel) including spin- 
orbit coupling. The energy is given relative to the valence- 
band maximum E v . As apparent in the case without SOC, 
mainly in-plane states contribute to the valence-band maxi¬ 
mum near K (Mo d x 2_ y 2 and d xy , Se p x / y ), while the max¬ 
imum near F is formed by out-of-plane states (Mo d z 2 , Se 
p z ). This also holds in the SOC case where the valence-band 
states with mainly in-plane character can be found near K 
(j = 5/2, rrij = ± 5/2, ±3/2 and j = 3/2, rrij = ± 3/2) while 
those with more out-of-plane character can be found near F 
(rrij = ± 1/2 for both j = 5/2 and j = 3/2, see Ref. [Hoi for the 
states in terms of spherical harmonics). On the other hand, 
the conduction-band minimum near K has mainly contribu¬ 
tions from out-of-plane Mo and in-plane Se states (without 
SOC: Mo d z 2 and Se p x / y , with SOC: Mo rrij = ±1/2 for 
both j — 5/2 and j = 3/2). 


MoSe 2 



FIG. 2. Band structure for bi- and trilayer MoSe 2 . The energy 
is given relative to the valence-band maximum E v . Molyb¬ 
denum diselenide changes (as all TMDs investigated in this 
paper) from a direct-band-gap semiconductor to an indirect 
one when the number of layers is increased. 


gap r —> K. Also for other TMDs one can find different 
result s 46 ' 59 ! 61- — . For a comprehensive review of the the¬ 
oretical papers see also Ref. H3 and references therein. 
Molybdenum ditelluride is especially peculiar, since the 
calculations show that its valence-band maximum is lo¬ 
cated at K even in the trilayer case. Furthermore, we 
find that the difference between the maximum at F and 
K is only 26 meV in bulk MoTe 2 which is in agreement 
with the experimental results in Ref. [tbil . 

In section UlI Dl we will see that, e.g., the varying dif¬ 
ference between the conduction-band minimum at K and 
Q for the different TMDs will also lead to a different 
thickness of the conductive channel for electron doping 
in an FET setup, while nearly all TMDs will behave sim¬ 
ilarly under hole doping. However, before investigating 
the changes under field-effect doping we want to focus on 
another problem that can make it difficult to dope a 2D 
system - the quantum capacitance. 


cordingly, the maximum (minimum) at T (Q) shift up 
(down) in energy which eventually leads to an indirect 
band gap between T and Q (see also Ref. [49| for an in- 
depth analyses). The transition between the indirect gap 
r —> K and T —> Q happens at different number of layers 
and occurs either in the bilayer case (MoSe 2 ), in the tri¬ 
layer case (WS 2 and WSe 2 ), or in the bulk limit (M 0 S 2 , 
MoTe 2 ). Depending on the level of theory that was used, 
one can also find very different results in literature when 
this transition occurs. Most calculations were done for 
M 0 S 2 for which the transition either already occurs for 
the bilaye r 46 : 55 : 56 or with larger number of layers^ - —. In 
fact, Ramasubramaniam et al. have shown in Ref. [60l that 
using the experimental bulk distance between the layers 
in the bilayer case also leads to an indirect band gap T —>■ 
Q while relaxation using PBE±D2 leads to an indirect 


B. Quantum capacitance 

A prominent example in which the quantum ca¬ 
pacitance hinders the doping via field-effect setup is 
graphene. Due to its linear dispersion relation at the 
K points, the charge that can be induced in an FET 
setup is much smaller than the corresponding charge at 
the gate electrod o ^^. Thus, doping concentrations ex¬ 
ceeding 10 13 cm -2 are hardly achievable using common 
dieletrics such as SiC >2 or HfC> 2 . Similarly, the quan¬ 
tum capacitance in nanolayers of TMDs could also re¬ 
duce the amount of induced charge. In the following we 
want to show that TMDs are however quite different from 
graphene and that the quantum capacitance is not rele¬ 
vant in their case as soon as the Fermi energy is within 
the conduction or valence band. In fact, experimen¬ 
tally doping concentrations in the order of 10 14 cm -2 are 
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FIG. 3. (a) Schematic illustration of an FET setup in which 
the 2D metallic system is separated from the gate electrode 
by a dielectric with dielectric constant e ox of thickness d ox . 
(b) Equivalent circuit for the overall capacitance seen at the 
gate electrode. 


TABLE IV. Geometrical capacitances for parallel-plate ca¬ 
pacitors with 100 nm SiC >2 (relative permittivity e r = 3.9), 
10 nm HfC >2 (e r ~ 20), and an 1 nm thick ionic-liquid (IL) 
FET (e r « 15), quantum capacitance for electron doping of 
the K valley (m « 0.5mo) and hole doping of the T valley 
(m « m o) of M 0 S 2 (v = 4), and the overall capacitance of 
a setup as shown in Fig. [3|b). All capacitances are given in 
units of gF cm -2 . 

_= 133.9_Cg = 267.7 

Csio 2 = 0.035 C = 0.035 C = 0.035 

C H fo 2 = 1-771 C = 1.748 C = 1.759 

Gil = 13.28 C = 12.08 C = 12.65 


possible using ionic-liquid based FET a 10 ' 16 ' 41 ’ 67 ' 68 (c/., 

Tab. HIIll. 

The term “quantum capacitance” was first used by 
Serge LuryiS* in order to develop an equivalent circuit 
model to describe the incomplete screening of an electric 
field by a 2D electron gas. When a 2D metallic system 
is contacted by a gate electrode (separated by a dielec¬ 
tric as shown in Fig. [Sja)), the electric field generated 
by the charges on the dielectric surface leads to a shift 
of the Fermi level of the 2D metal. This effect results in 
a modified capacitance with respect to the geometrical 
capacitance C ox per unit surface 

Cqx - E-OX&0 d ox } ( 1 ) 


obtained for a capacitor having a dielectric constant e ox 
and thickness d ox . 

As shown in Fig. [3jb) , the quantum capacitance Cq is 
in series with that of the dielectric, namely 


1 _ 1 1 
Wox '-'Q 


( 2 ) 


Under the simplified assumption that both C ox and Cq 
are independent of the applied gate voltage V G , the 
charge induced in the 2D system uq can be written as: 


0 V Vg e (V v ,V c ) 

n Q = l C(V v - V g ) v V G < 14 (3) 

( C(V c -V G ) VV G >V C 

where V v and V c are the onset potentials to fill the 
valence-band minimum or conduction-band maximum, 
respectively. Thus, when Cq C ox -, we have 

C = C ox (i + ^) -Cox- (4) 

and we regain the classical result, i.e., the charge that 
can be induced in the 2D system depends only on the 
applied gate voltage and the capacitance between gate 
and sample. 

As shown in Refs. SzO, the quantum capacitance in 
this case and for gate voltages larger (smaller) than the 


onset potential of the conduction band (valence band) is 
given by 


Cq 


g s g v mg 2 
nh 2 


(5) 


where g s and g v , m, and q are the spin and valley degen¬ 
eracies, the effective mass, and the charge, respectively. 
Depending on (i) the dielectric thickness d ox , (ii) the ef¬ 
fective mass of the 2D metallic system, and (iii) the num¬ 
ber of valleys v = g s g Vl the quantum capacitance Cq can 
be relevant or not. 

Table Ml shows the total capacitance, geometrical ca¬ 
pacitance for typical gate dielectrics, and the quantum 
capacitance for electron doping of the K valley or hole 
doping of the T valley of M 0 S 2 (the effective masses of 
the different valleys were taken from Ref. l46l ). As the 
geometrical capacitance of a parallel-plate capacitor in¬ 
creases with decreasing thickness of the dielectric, an 
ionic-liquid (IL) FET with an 1 nm thick electric double 
layer (inner Helmholtz plane) shows the largest deviation 
of the total capacitance C from the geometrical capaci¬ 
tance Cil ^ even for a thin dielectric layer of lOnm HfC> 2 , 
C ~ 0.986 CHfo 2 - Yet, for TMDs the number of valleys 
increases with increasing doping (valleys at K and Q or 
K and T for electron or hole doping, respectively), which 
also leads to a considerably larger DOS at the Fermi en¬ 
ergy than for a single quadratic band. Thus, the quantum 
capacitance Cq is further increased^ which in turn leads 
to C ~ C ox . Similar results have been found by Nan Ma 
and Debdeep Jena in Ref. [Tl| who also provide a detailed 
description of the low-doping regime n < 10 13 cm -2 


C. Structural changes under field-effect doping 

The doping via FET setup has only a minor influence 
on the structure of the TMDs - the changes are much 
smaller than for, e.g., ZrNCl^. This is due to the weak 
polarity of the bond between the transition-metal and 
the chalcogen. The largest change can be found for the 
layer thickness (z component of the chalcogen-chalcogen 
vector) of the layer closest to the charged plane repre¬ 
senting the gate electrode - the layer thickness increases 
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FIG. 4. (left) Structure of 3 different structural phases of 
monolayer M 0 S 2 . The H polytype is the one found in the 
bulk compound where the coordination of the molybdenum 
(gray) is trigonal prismatic. The T polytype with octahe¬ 
dral coordination can change to the T’ polytype in which the 
molybdenum atoms form zig-zag chains. The sulfur atoms 
are shown in yellow, (right) With increasing electron doping 
the difference between the total energy of 1T/1T’ structure 
(^x-MoS 2 ) and the total energy of the 1H polytype (-Eih-MosJ 
decreases. For a doping larger than n = —0.44 e/unit cell 
(n « —4.96 • 10 14 cm -2 ) 1 T’-MoS 2 is the lowest-energy struc¬ 
ture. 

by ~ 0.06 A for a large electron doping of n = — 0.3e/unit 
cell (n« - 3.16- 10 14 cm -2 ) and decreases by ~ 0.02A for 
a large hole doping of n = +0.3 e/unit cell. This change 
is mainly due to the increase/decrease of the chalcogen- 
transition-metal bond length of those being closest to the 
gate - we find ~ +0.04 A for n = —0.3 e/unit cell and 
~ —0.02 A for n = +0.3 e/unit cell. Accordingly, there is 
also a small change in the angle between the first chalco- 
gen, the transition-metal, and the second chalcogen of up 
to +0.9° (—0.4°) for large electron (hole) doping. Note 
that even if the structural changes seem to be small it 
is still important to relax the system in the FET setup. 
Otherwise, the band structure can be quite different for 
doping |n| > 0.15 e/unit cell as exemplified in Fig. 24 in 
the appendix. 

Even if the internal structure changes only slightly un¬ 
der FET doping, electron doping could induce a phase 
transition where the structure of the full nanolayer sys¬ 
tem is altered. It is well known that lithiu m 72 ! 73 or 
potassium^ intercalated M 0 S 2 can undergo a phase 
transition in which the Mo coordination changes from 
a direct-band-gap, semiconducting, trigonal-prismatic 
structure (labeled 2H, “2” as there are two layers in the 
unit cell) to a metallic, octahedral one (IT). It has been 
found experiment all y 75 1 76 and shown theoreticall y 77 ' 78 , 
that also monolayer M 0 S 2 can undergo this phase transi¬ 
tion. In the calculations, for a high electron doping by H 
(Ref. O or Li (Ref. |78|) adsorption of n « —0.35 e/unit 
cell or n = —0.44 e/unit cell, respectively, the octahedral 
phases such as the IT or IT’ phases become more stable 
than the 1H phase. In the IT’ phase the molybdenum 
forms zig-zag chains like tungsten in WTea^r— . Such a 
transition was also found experimentally by rhenium dop¬ 


ing of WS 2 nanotubes^ and monolayer MoS^, and by 
transfer of hot electrons generated in gold nanoparticles 
to monolayer MoS^- In order to determine if the FET 
setup can lead to such a phase transition for electron dop¬ 
ing, we compare in Fig. [4] the total energy for the mono- 
layer structures of IT-M 0 S 2 and 1 T’-MoS2 (^x-MoS 2 ) 
with the total energy of the 1H polytype (•E'ih-MoS 2 )' We 
find that the IT’ structure becomes more stable for elec¬ 
tron doping larger than n = —0.44 e/unit cell in close 
agreement with the results of Refs. I77lf78l . Thus, it seems 
that the interaction between the H/Li atoms and the 
M 0 S 2 layer has only a minor influence on the phase sta¬ 
bility as the transition occurs in our FET setup at the 
same doping. We also calculated the energy difference be¬ 
tween lT’-WSe 2 and lH-WSe 2 , (^. W Se 2 - ^-wsj, 
for a few electron-doping concentrations and found that 
for a concentration of n = —0.35 e/unit cell the IT’ poly¬ 
type becomes more stable by 63 meV. In the following, 
we will thus only consider doping of the H polytype with 
electron concentrations n > —0.35 e/unit cell as it is the 
most stable structure found in nature and is often used to 
prepare the samples by the mechanical-cleavage method. 

D. Band structure in FET setup 

In the following section, we want to investigate the 
influence of field-effect doping on the electronic proper¬ 
ties of the TMDs. The doping via FET setup changes 
the band structure considerably as exemplified in Figs. [5] 
and [6l which show the band structures for different elec¬ 
tron and hole doping levels for mono- and trilayer M 0 S 2 , 
respectively. In the appendix we also demonstrate that it 
is important to correctly model the FET setup by com¬ 
paring the band structures of mono- and trilayer MoSe 2 
calculated with a compensating jellium background to 
those calculated with our method (Figs. 25, 26). Fur¬ 
thermore, we also provide the band structures for more 
doping levels and all the other TMDs in the appendix 
(in total more than 250 calculations, Figs. 28-63). We 
summarized the evolution of the band structure with in¬ 
creasing doping in the left panel of Figs. [7] [TT] which show 
the position of the different band extrema with respect to 
the Fermi energy. Additionally, the right panel in those 
figures shows the relative amount of doping charge per 
valley given by 

£ 1 k< £ 2 

na = jrY Y i^.ki 2 ■ (6) 

“ keOc i 

Here a = {F,K,Q}, f l a defines the subset of k points 
which are closer to, e.g., a = T than to any a A I\ 
N a is the total number of those k points, and k is the 
eigenenergy for band i at k. The interval [ei, £ 2 ] is always 
chosen such that the probability density is integrated be¬ 
tween an energy within the former band gap E m and the 
Fermi energy Ep of the doped system, i.e., [Ep, E m ] and 
[E m ,Ep] for hole and electron doping, respectively. 
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Monolayer M0S2 


Electron doping 


Hole doping 





FIG. 5. Band structure for different FET induced doping of monolayer M 0 S 2 . The figures in the left column are for electron 
doping while the right column shows the hole doping case. For monolayer M 0 S 2 mainly the valleys at K are filled and only for 
a high doping of n as ±1.69 • 10 14 cm -2 (n = ±0.15 e/unit cell) a small amount of charge is in the maximum at F (hole doping) 
or in the minimum at Q (electron doping). The band structures for more doping levels and all the other TMDs can be found 
in the appendix (Figs. 28-63). 


Electron doping 

For n-type doping of monolayer M 0 S 2 (as for all mono- 
layer TMDs) the doping charge first occupies the extrema 
at K. For small doping (as long as only one valley is 
doped) the bands are rigidly shifted. However, as soon as 
a second valley is close to the Fermi energy, the doping 
cannot be described by a rigid shift of the bands any¬ 
more. For electron doping, the down shift of the bands 
at K slows down and, as the valley at Q starts to get 
occupied, is eventually reversed into an up shift. Finally, 
for high electron doping, the K valley is unoccupied again 


and the doping charge is solely localized around Q. Com¬ 
paring our results of the changes in the conduction band 
for electron doping of monolayer M 0 S 2 (Figs. [H [7]) with 
literature shows that it is important to correctly model 
the system - while the authors of Ref. [H find an up shift 
of the Q valley with increasing electron doping, we see a 
down shift. The opposite shift in Ref. [H might be due 
to the free-electron states at T (i.e., the states in the 
vacuum between the repeated images, c/., also Ref. [89.) 
which approach the Fermi en ergy with increasing dop¬ 
ing. Also the authors of Refs. [20P23ll89l find a down shift 
of the Q valley further supporting our results even if in 


































Trilayer M0S2 


Electron doping Hole doping 





FIG. 6. Band structure for different FET induced doping of trilayer M 0 S 2 . The figures in the left column are for electron 
doping while the right column shows the hole doping case. In contrast to the monolayer case, the doping at F/Q (hole/electron 
doping) is more important in trilayer M 0 S 2 . The band structures for more doping levels and all the other TMDs can be found 
in the appendix (Figs. 28-63). 


those works the asymmetric electric held in an FET has 
not been taken into account. The amount of doped elec¬ 
trons needed to have the charge completely localized at Q 
depends on the TMD ( i.e ., the initial energy difference 
between the minimum at K and Q) and is larger than 
n = — 2.2-10 14 cm -2 (n ss —0.2e/unit cell): the transition 
occurs for M 0 S 2 at n ~ —3.83 • 10 14 cm -2 , for MoSe 2 at 
n « —2.5-10 14 cm -2 , for MoTe 2 at n « — 2.22-10 14 cm -2 , 
for WS 2 at n « —3.41 • 10 14 cm -2 , and for WSe 2 at 
n « —3.31 • 10 14 cm -2 . Please note, that using LDA 
could slightly change these results as it leads to a smaller 
unit cell. The compressive in-plain strain would then re¬ 
duce the difference between the conduction-band minima 


at K and 045,46,48 

and thus decrease the doping-charge 
concentration needed to solely dope the valley at Q (c/., 
Fig. 27). 

In multilayer M 0 S 2 , WS 2 , and WSe 2 first the valley at 
K is doped and both valleys at K and Q are occupied until 
n « —3.3 • 10 14 cm -2 while in bi- and trilayer MoSe 2 and 
MoTe 2 the order is reversed: electrons first occupy the 
valley at Q and the doping at K is always smaller. This 
is due to the minimum at Q being lower in energy than 
the one at K in the undoped system (see Figs. 0 [9]). For 
a doping of n < —2.1 • 10 14 cm -2 for MoSe 2 the K valley 
is even unoccupied. Yet, one can expect that for small 
electron doping (|n| < 10 13 cm -2 ) of thick samples (more 





























than 3 layers) of M 0 S 2 , WS 2 , and WSe 2 the electrons will 
also first occupy the Q valley as this valley is lowered in 
energy with increasing number of layers. 
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Hole doping 

For p-type doping of the monolayer TMDs, the doping 
charge first occupies the extrema at K. However, in con¬ 
trast to the electron-doping case, in the high-hole-doping 
limit (n > +0.2 e/unit cell, n > +2.1 • 10 14 cm -2 ) of 
the monolayer TMDs both valleys at T and K are oc¬ 
cupied - the relative amount of doping charge in the T 
valley is even higher than that in K. The transition when 
the r valley is more occupied than the K valley again 
depends on the TMD, i.e., on the initial energy differ¬ 
ence between the valence-band maxima. In the case of 
M 0 S 2 for example even a small doping is enough to dope 
more the T valley. The only exception from this picture 
is MoTe 2 for which K is always more occupied even for 
n = +0.35e/unit cell (n ~ +3.23- 10 14 cm -2 , c/., Fig. [9]). 

Hole doping of the bilayer and trilayer TMDs is again 
very similar: in all investigated compounds (except 
MoTe 2 ) first the conduction-band maximum at T is oc¬ 
cupied, while for higher doping also K starts to get filled. 
Most interestingly, in some cases the second band at K 
is never occupied. For trilayer WS 2 and WSe 2 it is even 
pushed down in energy, effectively increasing the split¬ 
ting of the spin-orbit-split bands. For WSe 2 it is possible 
to achieve even higher hole-doping concentrations than 
shown in Fig. HU by using an IL-FET—. We thus calcu¬ 
lated for trilayer WSe 2 also some higher doping cases (up 
to n = +1 e/unit cell, or n ~ +1 ■ 10 15 cm -2 ) in order 
to understand what might happen for such a high dop¬ 
ing. Our calculated band structure is very different from 
the one in Ref. [TH which was calculated without proper 
treatment of the FET setup. For a doping of one hole per 
unit cell the second band at K is pushed down below the 
Fermi energy and also the first band at K is lowered. The 
band structure in Fig. [12] shows that the former band gap 
is closed and the first two bands are nearly unoccupied. 
However, since these are bands localized on the first layer 
close to the gate, we believe that for such a high doping 
the ions of the ionic liquid might start to interact with 
WSe 2 - In our simplified model without inclusion of the 
full dielectric it is difficult to prove this statement and we 
will thus leave this interesting problem for future investi¬ 
gations and will concentrate here on lower doping values. 


Conductive channel 

Another important property in order to understand 
different experiments on FET doping on TMDs is the 
doping-charge layer thickness, i.e., the size and shape 
of the conductive channel created and influenced by the 
gate voltage. To visualize the conductive channel we cal¬ 


culated the planar-averaged doping-charge distribution 
along z, 



e 

^2 D Nk 


£l<Ei,k<£2 

dA IV+k(r)| 2 (7) 

i,k 


Here is the unit cell area and Nk is the total number 
of k points. The interval [ei+ 2 ] is defined as above, i.e., 
[Ep,E m ] and \E m ,Ep\ for hole and electron doping, re¬ 
spectively. This property not only reveals the thickness 
of the conductive channel but also the relative distribu¬ 
tion among the different layers. 

In our calculations for low hole doping the charge is 
nearly evenly distributed between the first two layers 
with only small contributions at the third layer as can be 
seen in the upper panel of Fig. [13] Increasing the doping 
(i.e., increasing the gate voltage, lower panel of Fig. m 
the charge is more and more localized on the layer closest 
to the gate with a negligible amount of holes on the third 
layer. In the case of electron doping of multilayer M 0 S 2 , 
the charge is localized on the first layer as only the K val¬ 
ley is doped (small contribution of the out-of-plane states 
of sulfur). Only for higher doping when the minimum at 
Q is occupied a small amount of charge can be found on 
the second layer (c/., lower panel of Fig. [13] and right 
panel of Fig. 0. One can also see that the asymmetry 
in the doping-charge distribution is more pronounced for 
hole doping and that the system in this case moves closer 
to the dielectric. Furthermore, for low hole doping the 
doping charge on the second layer is even slightly larger 
than the one on the first. 

Figure [14] summarizes the doping-charge distribution 
for all bi- and trilayer TMDs. In the hole-doping case all 
TMDs behave similarly except MoTe 2 : for low doping the 
holes are delocalized over the first two layers with only 
small contributions in the third layer and, thus, the bi- 
and trilayer systems are nearly the same. Increasing the 
doping leads to stronger localization of the charge within 
the first layer and an effective narrowing of the conduc¬ 
tive channel. One can also easily understand why MoTe 2 
behaves differently: in all multilayer TMDs first the val¬ 
ley at r is doped, since it is the valence-band maximum, 
while in MoTe 2 the K valley is the maximum. As we 
have seen in section ITlI Al the states close to T have large 
out-of-plane contributions of both the transition metal 
and the chalcogen atom, while the valley at K is com¬ 
posed only of in-plane states. Increasing the hole doping 
however also leads for MoTe 2 to a small doping of the 
r valley. Thus, the amount of charge within the second 
layer increases slightly in the beginning. The small kink 
for trilayer MoTe 2 close to n = +0.46 ■ 10 14 cm~ 2 is due 
to both bands at T and K being close to the Fermi en¬ 
ergy. In this low doping limit (per valley) the calculation 
would require an infinite number of k points to fully con¬ 
verge the results. Further increase of the doping leads 
as for all TMDs to a larger screening of the electric field 
and therefore to the stronger localization within the first 
layer as can also be seen in the lower panel Fig. [13] 
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FIG. 7. Position of the different band minima/maxima Ei (left panel) with respect to the Fermi level and relative amount of 
doping charge per valley n a (right panel) as a function of doping for mono-, bi-, and trilayer M 0 S 2 . “Q” labels the conduction- 
band minimum half-way between K and T. In each graph, the scale and the units for the lower x axis are given in the lowest 
graph, while those of the upper x axis are given in the uppermost graph. Two different line styles for the two spin-orbit-split 
conduct ion-band minima at K were used in order to enhance the readability. Lines are guides for the eye. 


For electron doping we can divide the different TMDs 
into two different classes: (i) those in which the con¬ 
ductive channel for low doping [n ~ — 10 13 cm~ 2 ) has a 
thickness of just one layer (M 0 S 2 , WS 2 , WSe 2 ) and (ii) 
those with a three-layer-thick channel (MoSe 2 , MoTe 2 ). 
Using the results of Figs. [7][TT] one can see that in the 
TMDs of class (i) initially the K valley is doped while 
in (ii) the Q valley is occupied. Since the chalcogen 
states close to the conduction-band minimum at K have 
mainly in-plane character (in contrast to the transition- 
metal states which have d z 2 character), the hybridization 
between the layers is small and the electrons are more lo¬ 
calized within the first layer. The chalcogen states close 
to Q on the other hand have a large out-of-plane contri¬ 
bution which leads to a stronger hybridization between 
the layers. With increasing doping the electric field of the 
gate is more and more screened and the size of the con¬ 
ductive channel reduces to one layer. Furthermore, one 
can also understand why the tungsten dichalcogenides 
have a steep increase of the channel thickness in the be¬ 
ginning while this is not the case for M 0 S 2 : first, the 
difference between the conduction-band minimum at K 
and Q is much smaller in multilayer WS 2 /WSe 2 than in 
M 0 S 2 and second, a small electron doping can also re¬ 
sults in an effective separation of the single (doped) layer 
from the multilayer system. The difference between the 
conduction-band minimum at K and Q for monolayer 


M 0 S 2 is however twice as large as in the tungsten systems 
(~ 300 meV for MoS 2 and ~ 150 meV for WS 2 /WSe 2 ). 
Thus, the valley at the Q point is doped much earlier in 
WS 2 /WSe 2 than in M 0 S 2 . 

Number of occupied bands 

It is important to note, that the thickness of the 
doping-charge distribution, the number of occupied 
bands at a given k point, and the number of TMD layers 
(Le., the system size) are uncorrelated. Indeed, as can 
be seen in Fig. [15l the total spin-valley degeneracy v can 
be quite similar for different number of layers, whereas 
the doping charge is localized on one or two layers as 
seen in the previous section. Here the total spin-valley 
degeneracy has been calculated by counting the number 
of valleys within the interval [ei,£ 2 ] as defined above: 

= =^29s9v- ( 8 ) 

a. 

Here gf and < 7 “ are the spin and valley degeneracies of the 
valley at k point a = {F,K, Q}. For electron doping v 
is much higher than in the hole doping case as the valley 
degeneracy for the conduction-band minimum close to Q 
is g v = 6 (spin degeneracy g s = 1). Thus, as soon as 
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FIG. 8. Position of the different band minima/maxima Ei (left panel) with respect to the Fermi level and relative amount of 
doping charge per valley n a (right panel) as a function of doping for mono-, bi-, and trilayer MoSe 2 - “Q” labels the conduction- 
band minimum half-way between K and T. In each graph, the scale and the units for the lower x axis are given in the lowest 
graph, while those of the upper x axis are given in the uppermost graph. Lines are guides for the eye. 
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FIG. 9. Position of the different band minima/maxima Ei (left panel) with respect to the Fermi level and relative amount of 
doping charge per valley n a (right panel) as a function of doping for mono-, bi-, and trilayer MoTe 2 - “Q” labels the conduction- 
band minimum half-way between K and T. In each graph, the scale and the units for the lower x axis are given in the lowest 
graph, while those of the upper x axis are given in the uppermost graph. Lines are guides for the eye. 
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FIG. 10. Position of the different band minima/maxima Ei (left panel) with respect to the Fermi level and relative amount of 
doping charge per valley n a (right panel) as a function of doping for mono-, bi-, and trilayer WS 2 . “Q” labels the conduction- 
band minimum half-way between K and T. In each graph, the scale and the units for the lower x axis are given in the lowest 
graph, while those of the upper x axis are given in the uppermost graph. Lines are guides for the eye. 
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FIG. 11. Position of the different band minima/maxima Ei (left panel) with respect to the Fermi level and relative amount of 
doping charge per valley n a (right panel) as a function of doping for mono-, bi-, and trilayer WSe 2 - “Q” labels the conduction- 
band minimum half-way between K and T. In each graph, the scale and the units for the lower x axis are given in the lowest 
graph, while those of the upper x axis are given in the uppermost graph. Lines are guides for the eye. 
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Trilayer WSe 2 



FIG. 12. Band structure of trilayer WSe 2 for a high hole 
doping of n « +1.03 • 10 15 cm~ 2 (n = +1 e/unit cell). 



FIG. 13. Planar-averaged doping-charge distribution along 
^ for trilayer M 0 S 2 for a doping of n « ±2.25 ■ 10 13 cm -2 
(n = ±0.02e/unit cell, upper panel) and n ps ±1.69T0 14 cm -2 
(n = ±0.15 e/unit cell, lower panel). The dashed line within 
the graph indicates the end of the barrier potential, while the 
sketch above shows the position of the atoms (gray - Mo, 
yellow - S). 


the valley at Q is doped, v increases drastically by 6 or 
12 depending on whether only one or both spin-orbit- 
split bands are filled. For some TMDs a high electron 
doping of n < —3.5 • 10 14 cm -2 leads to a large lowering 
of the minimum at Q - so much that it is actually not 
a valley anymore but a ring around T. This can be seen 
by the minimum between T and M which appears in the 
band structure (c/., Figs. [5J [G] and band structures in 
the appendix). In this case, we do not count it as six 


Electron doping Hole doping 



n [ I0 14 cm" 2 ] 


FIG. 14. Relative doping charge per layer for bi- and trilayer 
of the different dichalcogenides with increasing FET doping. 


independent valleys but as one. 

In contrast, for high hole doping (n > +2-10 14 cnG 2 ) of 
multilayer TMDs the number of occupied valleys is either 
v = 6 or v = 8 as only bands at K ( g v = 2, g s = 1) and 
r ( g v = 1, g s = 2) are doped. In the monolayer systems 
often two valleys are doped — the valence-band maximum 
at K and either the second spin-orbit-split band at K or 
the valence-band maximum at T. Accordingly, the dif¬ 
ference between the valence-band maximum at K and 
at r determines the doping-charge concentration needed 
in order to dope two valleys at different points in the 
BZ. The spin-degenerate maximum at T is occupied for 
monolayer M 0 S 2 , MoSe 2 , WS 2 , and WSe 2 for hole doping 




















































14 


Electron doping 
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FIG. 15. Total spin-valley degeneracy v as function of dop¬ 
ing for the mono-, bi-, and trilayer system of all investigated 
TMDs. The total spin-valley degeneracy has been calculated 
by summing the spin and valley degeneracies of all doped val¬ 
leys. Lines are guides for the eye. 


larger than +1 • 10 13 cm -2 , +5-10 13 cm -2 , +2- 10 13 cm -2 , 
and + 8 - 10 13 cm -2 , respectively. Again monolayer MoTe 2 
is exceptional due to the large difference of w 600meV be¬ 
tween the maximum at K and T. Just for a very high hole 
doping of n ~ +3.23 • 10 14 cm -2 (n = +0.35 e/unit cell) 
the second band at K is occupied and a “ring” around 
and close to T appears (which is again counted as 1 ). 

Up to now we have focused on the changes in the elec¬ 
tronic structure of the different TMDs with increasing 
doping and we saw that for, e.g., high electron doping 


the charge is mainly localized around the Q point. We 
now want to investigate how the amount of doping charge 
and thus the number of occupied bands and the thickness 
of the conductive channel is determined experimentally. 


E. Hall-effect measurements 

In order to determine the doping charge in a sample, 
one commonly performs a Hall experiment as the inverse 
Hall coefficient is directly proportional to the charge- 
carrier density n in the case of parabolic, isotropic bands. 
This, however, is already a crude approximation which, 
as we will see below, can lead to large differences be¬ 
tween the charge-carrier density thus calculated and the 
real density within the sample. 

We closely follow the work of Madsen and Singh^ 
and sketch the calculation of the Hall coefficient (or 
more specifically, the Hall tensor Rijk ) within Boltzmann 
transport theory. In the presence of an electric field E 
and a magnetic field B, the electric current j can be writ¬ 
ten as 


ja — kT a gE/^ + &a(3~/R/3 By + • • * , (9) 

with the conductivity tensors cr a p and cr a( g 7 , and the in¬ 
dices denoting the spatial dimensions. Here and hence¬ 
forth, we will always adopt Einstein’s sum convention, 
according to which whenever an index occurs twice in 
a single-term expression, the summation is carried out 
over all possible values of this index. The Hall tensor is 
defined as 


Rijk (T]Ef) 


jpind 


• appl t -yappl 

J; 


( CT 1 )cj ( a X )i& ’ 

( 10 ) 


where is the electric field along direction j which is 
induced by the applied magnetic field B)( ppi and current 
jappZ a j on g direction k and i, respectively. 

Within the relaxation-time approximation the conduc¬ 
tivity tensors cr a p and cr a p^ for temperature T and chem¬ 
ical potential E F are given in the 2D case by^ 


<T a p ( T-E f ) = -^-2 f T i>k vv)'f 
x I" dfpp (T; £+ k ) 


de 


d 2 k, 


0+87 (T;Ep) = 


Q 


i.k 


(2tt) 


T i,k €r yuv v d v d 


( 11 ) 


dfEp [T] £*, k ) 


de 


d 2 k. 


( 12 ) 


Here, q = ±e is the charge of the charge carriers in band 
Ei k with momentum k, f FF (T;e) is the Fermi function 
f E p (T ; s) = (exp[(e-E F )/(k B T)] + l)- 1 , and e 7UI/ is the 
Levi-Civita symbol. The relaxation time 7 y k in principle 
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is dependent on both the band index i and the k vector 
direction. Furthermore, is the group velocity 


v 


i,k 

OL 


1 dsi,k 

h dk a 


and 



the inverse mass tensor 



J_ d 2 £j^ 
Ti 2 dkpdku 


(13) 


(14) 


To show that the inverse Hall coefficient R X yz is 
proportional to n, we start by assuming bands with 
quadratic dispersion. The dispersion relation for a 
quadratic, isotropic band in 2D is given by 


&i,k — 


n 2 k 2 

2 rrii ’ 


(15) 


with k 2 = k 2 + k 2 . The group velocity is then v^ ka = 
hk a /m,i while the mass tensor in Eq. m is for each band 
a diagonal matrix with to ^ 1 on the diagonal. In the zero 
temperature limit and assuming an i- and k-independent 
relaxation time r,.k = r[Ep) the conductivity distribu¬ 
tions in Eqs. CD and m are given by 

a 2 t 

<j aa (0;E F ) = V]- m, (16) 

^ rm 

l 

q 3 r 2 

&cx/3 7 ( 0 ? Bjr') — ^ ^ 2 ( 17 ) 

ttlz 


where rq is the charge-carrier density in band i. Finally, 
assuming that the magnetic field is applied perpendicular 
to the 2D system along z, we get for the Hall coefficient 

R xyz (0; E F ) = ^ 2 ■ (18) 

9 ( z2i m i n i) 

Thus, only for valley-independent effective mass rm = m 
the inverse Hall coefficient is directly proportional to the 
doping-charge concentration n = n q. The results in 
Eqs. (USD, (1171) . and (USD also hold for a 3D system, how¬ 
ever with the important difference that the conductiv¬ 
ity in two dimensions is independent of the mass of the 
charge carriers, since the density n is proportional to to. 
Furthermore, the Hall coefficient is inversely proportional 
to the mass to in the 2D case. 

In the constant-relaxation-time approximation (r,;.k = 
t(Ef)) and for hexagonal symmetry (such as in the 
TMDs), the Hall tensor in Eq. (fToTl has only two inde¬ 
pendent coefficients^ (in-plane and out-of-plane compo¬ 
nent). However, it is important to remember that this 
simple equation for the Hall coefficient R xyz , Eq. (usd. is 
only valid as long as there are only bands with isotropic, 
quadratic dispersion and if 7 yk = t(Ef). For small 
doping this might be a good approximation but in the 
high doping case this approximation can break down - 
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FIG. 16. Ratio of the inverse Hall coefficient R X y Z to the 
doping charge n as a function of doping for the mono- (black, 
solid), bilayer (red, dashed), and trilayer (green, dash-dotted) 
of all investigated TMDs for a temperature of T = 300 K. 
Note also the different range of the ordinate in the case of 
MoTe 2 • 


especially, as soon as a minimum/maximum with non¬ 
quadratic dispersion starts to get filled. 

Figure Uni shows the ratio of the inverse Hall coefficient 
R~y Z (T- Ep) (calculated with the BoltzTraP code^ us¬ 
ing Eqs. (fill (fl^l) and assuming 7 y.k = r(Ep )) to the 
doping-charge concentration n, as a function of doping 
for the mono-, bi-, and trilayer systems for T = 300 K. 
The comparison between T = 0 K and T = 300 K which 
can be found in the appendix, Figs. 64-66, shows that 
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Monolayer MoTe 2 



FIG. 17. Hall coefficient R xyz as a function of doping for 
monolayer MoTe 2 and temperatures T = 300 K, T = 100 K, 
and T = 50 K. The inset shows the band structure for a 
critical doping of n « 1.48 • 10 14 cm -2 (n = +0.16e/unit cell). 
The specific band structure of MoTe 2 with the nearly linear 
dispersion along K —»• T and the changing sign of the effective 
mass when increasing the distance to K leads to the large 
difference between R xyz and n. 


the temperature has only a minor influence on the in¬ 
verse Hall coefficient except for MoTe 2 (which we will 
also discuss separately). 

For most of the studied systems, the inverse Hall co¬ 
efficient R xyz (300 K; Ep) shows a strong deviation from 
the doping charge n for large doping - the doping-charge 
concentration calculated using the inverse Hall coefficient 
can be 1.5 times larger than the real concentration. In 
fact, R~y Z (300K; Ep) « n is true only in two cases (c/., 
Eq. ([18]) ): either (i) if mainly one valley with parabolic 
dispersion is doped or (ii) if the doping charge is split 
between several valleys with quadratic bands which how¬ 
ever have similar effective masses. 

Case (i) holds only for small hole doping (n < +5 • 
10 13 cm -2 ) of all multilayer TMDs as mainly the V val¬ 
ley is doped. The deviation of the inverse Hall coeffi¬ 
cient from the doping charge n increases with increasing 
doping of the valence-band maximum at K. Accordingly, 
R xyz (300 K; Ep) « n for a larger range of hole doping of 
multilayer M 0 S 2 than in the other TMDs as the doping at 
K is negligible. For high hole doping (n > +2-10 14 cm -2 ) 
the deviation is due to the non-parabolicity of the bands 
near K. If however the charge is split between T and K 
the difference can be explained by the much larger effec¬ 
tive mass of the T valley (i.e., in variance with case (ii)). 
Assuming np ss nx = n/2, the ratio of the inverse Hall 
coefficient to the doping charge (c/., Eq. (ITS]) ) simplifies 
to 


Rxyz (0; Ep) 1 


1 

2 


mr rriK 

9 ! 2 ~ ’ 

TOf. + m £■ 


(19) 


which is always smaller than 1 for mr ^ m k- This also 
explains why the ratio approaches ss 0.7 for small hole 
doping of monolayer M 0 S 2 . At T = 300 K n is split 
between both maxima which have however very different 
masses^ - mp = 3.524 mo and m k = 0.637 mo- 

The agreement between i?“^ z (300K; Ep) and n is much 
better for electron doping up to n « —2 ■ 10 14 cm -2 even 
if both conduction-band minima at K and Q are occu¬ 
pied. This is due to the similar effective masses of those 
two valleys. As in the hole-doping case the difference 
between i?Tyz(300 K; Ep) and n increase for larger dop¬ 
ing (n )$ —2 • 10 14 cm -2 ) which is due to the increasing 
non-parabolicity of the bands. 

The only case where the model of a 2D electron gas 
with constant relaxation time gives reasonable results for 
the doping-charge concentration within a large range of 
both electron and hole concentrations is multilayer M 0 S 2 . 
Since the agreement between R xyz (T ; Ep) and n can be 
much better in other systems as exemplified for CoSb 3 
in Ref. HD, the deviations shown in Fig. [TO] point out 
problems if the specific band structure is not taken into 
account. Once again, MoTe 2 is particularly interesting 
because the inverse Hall coefficient R xyz (300 K; Ep) is 
nearly three times bigger than n for a hole doping of 
n ~ +1.4 ■ 10 14 cm~ 2 . 

In order to understand the origin of this behavior of 
Rxyz I n f° r hole doping of monolayer MoTe 2 , we plot 
Rxyz as a function of temperature in Fig. 1171 As the 
temperature is reduced, the Hall coefficient R xyz {T; Ep) 
decreases and for T < 50 K it even changes the sign. The 
result for T = 50 K in Fig. [T7] indicates that in the range 
of +1.25• 10 14 cm -2 < n < +1.5 • 10 14 cm -2 the Hall coef¬ 
ficient R xyz changes twice the sign. The band structure 
for a doping of n s=s +1.48 • 10 14 cm -2 (inset in Fig. [T71) 
shows that the valence band has at least two inflection 
points close to the K point which cross the Fermi energy 
with increasing doping. Accordingly, the effective mass 
changes the sign and thus also the conductivity tensor 
(7 a fj-y and Rijk ■ Furthermore, the nearly linear dispersion 
along K —»• r leads tom + 0 and thus to the large differ¬ 
ence between R xyz and n. For bi- and trilayer MoTe 2 the 
deviations are smaller and at higher doping values. This 
is due to the finite contribution of the doping charge at 
the r point which leads to a smaller doping around K. 


F. Conductivity and DOS at the Fermi energy 

In the end, we want to briefly analyze the DOS at the 
Fermi energy Ep and the in-plane conductivity o+ x /t dr 
Fig.HH Both were calculated using the fitted band struc¬ 
ture of BoltzTraP. The DOS at the Fermi energy Ep in 
the left-hand panel shows, that the doping charge cannot 
always be described with quadratic, isotropic bands in 2D 
(as also shown above by the behavior of R xyz )- In this 
case, DOS (Ep) would be constant and would have steps 
as soon as another band crosses Ep. In fact, DOS (Ep) 
has steps and those can be related to crossing bands, but 


n 
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FIG. 18. DOS at the Fermi energy Ef (left panel, T = 0 K) and in-plane conductivity ff„/r (right panel, T = 300 K) as 
a function of doping for the mono- (black, solid), bilayer (red, dashed), and trilayer (green, dash-dotted) of all investigated 
TMDs. Note also the different range of the ordinate for DOS(-Ef) in the case of WS 2 and WSe 2 . Lines are guides for the eye. 


for high doping it can deviate from a simple 2D electron 
gas. 

For hole doping of monolayer M 0 S 2 , MoSe 2 , WS 2 , and 
WSe 2 the DOS at the Fermi energy is nearly constant 
as soon as the T valley is doped, i.e., for doping larger 
than n « +0.1 • 10 14 cm -2 , n « +1.5 • 10 14 cm -2 , n w 
+0.5 • 10 14 cm -2 , and n « +0.9 • 10 14 cm -2 , respectively. 
Also for hole-doping smaller than n < +1.1 • 10 14 cm -2 of 
multilayer M 0 S 2 the DOS is constant as only the T val¬ 
ley is doped in this regime. The non-constant behavior 
for hole doping of the other multilayer systems increases 
with increasing doping of the K valley: it is more pro¬ 


nounced for MoSe 2 and WSe 2 than for M 0 S 2 and WS 2 
(c/., Figs. ITl fTTI) . In the case of MoTe 2 , where mainly 
the valence bands at K are doped, a description with 2D, 
quadratic, isotropic bands completely fails. 

For n-type doping of all TMDs the DOS at the Fermi 
energy shows a quasi 2D behavior. It has steps when 
the conduction-band minimum at K or Q enters the bias 
window and is nearly constant in between. However, 
for larger electron doping the non-constant behavior in¬ 
creases. This is due to the stronger deviation of the spin- 
orbit-split conduction band at Q from a quadratic dis¬ 
persion. 
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The in-plane conductivity <t xx /t is another measure 
for the deviation from quadratic, isotropic bands in 2D. 
As can be seen in Eq. (fTfil) . for a perfect 2D electron 
gas the conductivity would be an increasing linear func¬ 
tion of the doping-charge concentration n. Most interest¬ 
ingly, the right-hand panel of Fig. fl8l shows that (Jxx/t 
weakly differs from one TMD to the other and from the 
monolayer to the multilayer case. Only for the tungsten 
dichalcogenides the conductivity has a small nonlinear 
component for n-type doping. This might be due to the 
stronger SOC of tungsten which leads also to larger de¬ 
viation of the Q valley from a quadratic dispersion. 

IV. CONCLUSIONS 

In this work, we have calculated from ab initio 
the structural, electronic, and transport properties for 
mono-, bi-, and trilayer TMDs in field-effect configura¬ 
tion. We have first investigated the structural changes 
of the TMDs under field-effect doping. We found that 
the internal structure is only slightly changed but that it 
is nevertheless important to fully relax the system. We 
also showed that high electron doping can induce a phase 
transition from the 1H to the IT’ structure. In accor¬ 
dance with literature this transition can however only 
occur for electron doping larger than n < —0.35 e/unit 
cell. Therefore, we concentrated on smaller doping of 
the H polytype as it is the most stable structure found 
in nature. 

The band structure and thus also the transport prop¬ 
erties can be changed considerably under field-effect dop¬ 
ing. We have shown that most TMDs behave similarly 
under hole doping while for electron doping they can be 
divided into two different classes: one in which the con¬ 
ductive channel has a width of approximately two layers 
for small doping (MoSe 2 and MoTe 2 ) and one in which 
the charge is localized within the first layer. This can 
be attributed to the relative position of the conduction- 
band minimum at K and Q in the multilayer TMDs. 
In the former class Q is lower than K and the doping 
charge first occupies the states in the Q valley. Since 
these states have a large p z contribution of chalcogen 
states, the hybridization between the layers is larger and 
the electrons are more delocalized. Additionally, for the 
tungsten dichalcogenides the difference between K and Q 
is smaller than for M 0 S 2 and thus, as the electron doping 
is increased, the charge rapidly starts to occupy states at 
Q. This leads to an increase of the width of the conductive 
channel to approximately two layers. For high electron 
doping only the Q valley is occupied in all investigated 
TMDs (also in the monolayer systems) and the width of 
the conductive channel is reduced to one layer. 

Under hole doping most TMDs behave similarly: in 


the monolayer case first the valence-band maximum at 
K is occupied while in the multilayer case it is the T 
valley. This also leads to the delocalization of the doping 
charge over more layers as the states at T have large 
chalcogen p z character. For large hole doping also in the 
monolayer case the doping at T is larger than those at 
K. Only exception from this picture is MoTe 2 in which 
the valence-band maximum for the undoped compound 
is always at K even in the multilayer case. Accordingly, 
the doping charge within the K valley is always larger 
than those at T. 

However, even if the thickness of the doping-charge 
distribution is approximately two layers for multilayer 
TMDs, the number of doped valleys can be comparable 
to the monolayer case. The main difference is between 
n-type and p-type doping: as the valley degeneracy for 
the conduction-band minimum at Q is g v = 6 , the total 
number of doped valleys can be as large as v = 16 while 
for hole doping the maximum is v = 8 . 

In the next part, we have seen that a Hall-effect mea¬ 
surement can often not directly be used in order to de¬ 
termine the charge-carrier concentration under the as¬ 
sumption of quadratic bands - the charge thus deter¬ 
mined can be up to 1.5 times larger than the real doping 
charge within the sample. For MoTe 2 the Hall coefficient 
Rxyz even changes the sign due to the changing curva¬ 
ture of the valence band. Thus, an interpretation based 
on parabolic bands would lead to an incorrect sign of 
the charge of the carriers. Even if this can only be seen 
at low temperatures, the Hall-effect measurement still 
largely overestimates the doping charge concentration if 
one does not take into account the specific band struc¬ 
ture. Only in the case of multilayer M 0 S 2 the inverse Hall 
coefficient is directly proportional to the doping-charge 
concentration for a large range of electron and hole dop¬ 
ing. 

In this work, we have shown that it is important to cor¬ 
rectly model an FET setup. The changes in the electronic 
and transport properties cannot be described with both 
the rigid doping and the uniform-background doping ap¬ 
proach. We provide not only a full database of electronic 
structure of mono- , bi-, and trilayer dichalcogenides as 
a function of doping, but also a mapping between the 
doping charge and the Hall coefficient. 
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Appendix A: Band structure of the undoped TMDs 



FIG. 19. Band structure of mono-, bi-, trilayer, and bulk M 0 S 2 . The arrow indicates the lowest-energy transition. 
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FIG. 20. Band structure of mono-, bi-, trilayer, and bulk MoSe 2 - The arrow indicates the lowest-energy transition. 
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FIG. 21. Band structure of mono-, bi-, trilayer, and bulk MoTe 2 - The arrow indicates the lowest-energy transition. 
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FIG. 22. Band structure of mono-, bi-, trilayer, and bulk WS 2 . The arrow indicates the lowest-energy transition. 
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FIG. 23. Band structure of mono-, bi-, trilayer, and bulk WSe 2 - The arrow indicates the lowest-energy transition. 















26 


Appendix B: Influence of structural relaxation in FET setup 


n = -0.35 e/unit cell n = -0.15 e/unit cell 





FIG. 24 . Comparison of the band structures of doped bilayer M0S2 with and without relaxation. The doping is indicated in the 
labels. For the fixed structures we used the relaxed geometry of undoped bilayer M0S2 and placed the system at approximately 
the same distance to the barrier potential as in the case of the fully relaxed system. 
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Appendix C: Influence of the asymmetric electric field in an FET 


Monolayer, n = -0.15 e/unit cell Monolayer, n = -0.02 e/unit cell 



Monolayer, n = +0.02 e/unit cell Monolayer, n = +0.15 e/unit cell 



FIG. 25. Comparison of the band structures of doped monolayer MoSe 2 calculated with a compensating jellium background 
(homogeneous-background doping, “h.b.d.”) and the proper FET setup. The doping is indicated in the labels. All systems 
have been fully relaxed. The difference is negligible for a small doping of n = ±0.02 e/unit cell, but for a higher doping of 
n = ±0.15 e/unit cell one can see clear differences not only in the relative amount of doping in the different valleys but also 
in the curvature of the bands, i.e., their effective masses. Furthermore, it is difficult to model the system properly with a 
compensating jellium background for electron doping larger than n = —0.15 e/unit cell due to the interlayer states crossing the 
Fermi energy. 
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FIG. 26. Comparison of the band structures of doped trilayer MoSe 2 calculated with a compensating jellium background 
(homogeneous background doping, “h.b.d.”) and the proper FET setup. The doping is indicated in the labels. All systems 
have been fully relaxed. The difference is again small for a low doping of n = ±0.02 e/unit cell. However, in contrast to the 
monolayer case, in the FET setup other states are close to the Fermi energy which will change the Hall coefficient for finite 
temperature. For a higher doping of n = ±0.15 e/unit cell one can see clear differences not only in the relative amount of 
doping in the different valleys but also in the curvature of the bands, i.e., their effective masses. The differences are also larger 
for high hole doping as shown in the lower right panel. 
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Appendix D: Influence of the exchange-correlation approximation 


Monolayer MoS 2 Monolayer MoS 2 , n = -0.10 e/unit cell 



Bilayer MoS,,, n = -0.10 e/unit cell Trilayer MoS,,, n = -0.10 e/unit cell 



FIG. 27. Comparison of the band structures of (doped) M 0 S 2 calculated using either PBE+D2 or LDA for the exchange- 
correlation energy. In the LDA case we used the final, fully-relaxed geometry of the PBE+D2 calculations except for the doped 
monolayer for which we also show the results for the LDA-relaxed unit cell (“LDA - fully relaxed”). One can see only minor 
differences in the relative occupation of the conduction-band minima at K or Q if the geometry of PBE+D2 is used. This is due 
to the different interlayer spacing in the case of PBE+D2 or LDA. Thus, as we used the geometry of the PBE+D2 calculations, 
LDA leads to minor changes in the bands which are doped. That this is only a result of the different interlayer spacing is 
further confirmed by the fact that the difference in the case of the monolayer is negligible. In the fully relaxed LDA case the 
differences are much larger which is due to the compressive strain as the LDA lattice parameter is 2.2% smaller than the one 
calculated using PBE+D2. 
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Appendix E: Band structure of the doped TMDs 


1. Molybdenum disulfide 
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FIG. 28. Band structure of monolayer M 0 S 2 for different doping as indicated in the labels. 
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n = +0.07 e/unit cell 



n = +0.08 e/unit cell 



FIG. 29. Band structure of monolayer M 0 S 2 for different doping as indicated in the labels. 
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n = +0.10 e/unit cell 


n = +0.15 e/unit cell 
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FIG. 30. Band structure of monolayer M 0 S 2 for different doping as indicated in the labels. 
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FIG. 31. Band structure of bilayer M 0 S 2 for different doping as indicated in the labels. 
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FIG. 32. Band structure of bilayer M 0 S 2 for different doping as indicated in the labels. 
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FIG. 33. Band structure of bilayer M 0 S 2 for different doping as indicated in the labels. 
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FIG. 34. Band structure of trilayer M 0 S 2 for different doping as indicated in the labels. 
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n = +0.02 e/unit cell n = +0.05 e/unit cell 



FIG. 35. Band structure of trilayer M 0 S 2 for different doping as indicated in the labels. 
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2. Molybdenum diselenide 
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n = -0.15 e/unit cell 
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n = -0.08 e/unit cell n = -0.05 e/unit cell 




n = -0.02 e/unit cell 



n = -0.01 e/unit cell 



FIG. 36. Band structure of monolayer MoSe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



FIG. 37. Band structure of monolayer MoSe 2 for different doping as indicated in the labels. 
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n = -0.08 e/unit cell n = -0.05 e/unit cell 




n = -0.02 e/unit cell n = -0.01 e/unit cell 




FIG. 38. Band structure of bilayer MoSe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell 



n = +0.35 e/unit cell 



FIG. 39. Band structure of bilayer MoSe 2 for different doping as indicated in the labels. 
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n = -0.35 e/unit cell 


n = -0.20 e/unit cell 



n = -0.02 e/unit cell n = -0.01 e/unit cell 




FIG. 40. Band structure of trilayer MoSe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell 



n = +0.35 e/unit cell 



FIG. 41. Band structure of trilayer MoSe 2 for different doping as indicated in the labels. 
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3. Molybdenum ditelluride 
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FIG. 42. Band structure of monolayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.12 e/unit cell n = +0.13 e/unit cell 




FIG. 43. Band structure of monolayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.135 e/unit cell n = +0.14 e/unit cell 



n = +0.17 e/unit cell 



n = +0.175 e/unit cell 



FIG. 44. Band structure of monolayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.18 e/unit cell 



n = +0.20 e/unit cell 



n = +0.35 e/unit cell 



FIG. 45. Band structure of monolayer MoTe 2 for different doping as indicated in the labels. 
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FIG. 46. Band structure of bilayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.10 e/unit cell 



n = +0.15 e/unit cell 



n = +0.20 e/unit cell 



n = +0.35 e/unit cell 



FIG. 47. Band structure of bilayer MoTe 2 for different doping as indicated in the labels. 
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FIG. 48. Band structure of trilayer MoTe 2 for different doping as indicated in the labels. 
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FIG. 49. Band structure of trilayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.04 e/unit cell n = +0.05 e/unit cell 



FIG. 50. Band structure of trilayer MoTe 2 for different doping as indicated in the labels. 
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n = +0.20 e/unit cell n = +0.35 e/unit cell 



FIG. 51. Band structure of trilayer MoTe 2 for different doping as indicated in the labels. 
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4. Tungsten disulfide 
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n = -0.05 e/unit cell 



n = -0.02 e/unit cell n = -0.01 e/unit cell 




FIG. 52. Band structure of monolayer WS 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell n = +0.35 e/unit cell 




FIG. 53. Band structure of monolayer WS 2 for different doping as indicated in the labels. 



































(E-E f ) [eV] (E-E p ) [eV] 


1.0 

0.8 

0.6 


n = -0.35 e/unit cell 


n = -0.20 e/unit cell 



n = -0.15 e/unit cell n = -0.10 e/unit cell 




n = -0.08 e/unit cell n = -0.05 e/unit cell 
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FIG. 54. Band structure of bilayer WS 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell 



n = +0.35 e/unit cell 



FIG. 55. Band structure of bilayer WS 2 for different doping as indicated in the labels. 
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FIG. 56. Band structure of trilayer WS 2 for different doping as indicated in the labels. 
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FIG. 57. Band structure of trilayer WS 2 for different doping as indicated in the labels. 
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5. Tungsten diselenide 
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n = -0.20 e/unit cell 




n = -0.15 e/unit cell 



n = -0.10 e/unit cell 



n = -0.08 e/unit cell n = -0.05 e/unit cell 




n = -0.02 e/unit cell n = -0.01 e/unit cell 




FIG. 58. Band structure of monolayer WSe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell 
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FIG. 59. Band structure of monolayer WSe 2 for different doping as indicated in the labels. 
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FIG. 60. Band structure of bilayer WSe 2 for different doping as indicated in the labels. 
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n = +0.01 e/unit cell n = +0.02 e/unit cell 



n = +0.20 e/unit cell n = +0.35 e/unit cell 




FIG. 61. Band structure of bilayer WSe 2 for different doping as indicated in the labels. 
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FIG. 62. Band structure of trilayer WSe 2 for different doping as indicated in the labels. 
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FIG. 63. Band structure of trilayer WSe 2 for different doping as indicated in the labels. 
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Appendix F: Hall measurements 


Electron doping Hole doping 



FIG. 64. Ratio of the inverse Hall coefficient R X y Z to the doping charge n as function of doping for all monolayer TMDs and 
temperatures T = 0 K and T = 300 K. 
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Electron doping Hole doping 



FIG. 65. Ratio of the inverse Hall coefficient R X y Z to the doping charge n as function of doping for all bilayer TMDs and 
temperatures T = 0 K and T = 300 K. 






























73 


Electron doping Hole doping 



FIG. 66. Ratio of the inverse Hall coefficient R X y Z to the doping charge n as function of doping for all trilayer TMDs and 
temperatures T = 0 K and T = 300 K. 




























